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This  report  deals  with  the  potential  boundary  layer  destabilizing 
action  of  machinery  generated  structure-borne  noise  in  submarine  hull 
plating.  Boundary  layer  transition  is  assumed  to  be  governed  by  the 
eigenvalues  of  the  Orr-Sommerf eld  eguation  and  the  structure-borne  noise 

field  is  taken  to  be  flexural  in  nature.  The  structure-borne  noise 2_ 
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'field  is  decomposed  into  propagating  wave  fields  and  near-fields  such  as 
would  surround  foundation  structures  or  other  hull  impedance  discontinu- 
ities. It  is  found  that  for  the  frequency  range  generally  associated 
with  structure-borne  noise,  say  30  - 10^  Hz,  the  propagating  wave  fields 
are  not  destabilizing  in  that  they  do  not  contain  components  which  are 
coincident  with  the  appropriate  eigenvalues  of  the  Orr-Sommerf eld  equation. 
By  contrast,  the  near-fields  will  contain  components  which  are  coincident 
with  the  potentially  destabilizing,  eigenvalues  of  the  Orr-Somnerf eld 
equation.  These  near-fields  are  interpreted  in  terms  of  equivalent 
propagating  fields  to  which  they  are  found  to  be  of  (relatively)  low 
level.  Available  experimental  results  are  discussed  in  terms  of  the 
submarine  environment  and  it  is  found  that  although  transition  may  be 
affected  by  structure-borne  noise  levels  less  than,  say,  the  expected 
free-stream  turbulence  levels,  the  measured  noise  levels  required  for 
destabilization  are  higher  than  one  expects  in  submarine  hull  plating. 

Finally,  a mechanical  analogue  to  the  Orr-Sommerf eld  equation  is 
presented  in  the  form  of  a vibrating  elastic  plate  resting  on  a locally 
reacting  foundation  of  specific  form.  In  the  analogue  Reynolds  number 
becomes  the  cube  of  the  slenderness  ratio  of  the  plate.  Some  aspects 
of  the  solutions  to  the  Orr-Sommerf eld  equation,  as  it  is  applied  to 
boundary  layer  transition,  are  interpreted  in  terms  of  the  analogue 
system. 
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SUMMARY 

This  report  deals  with  the  potential  boundary  layer  destabilizing 
action  of  machinery  generated  structure-borne  noise  in  submarine  hull 
plating.  Boundary  layer  transition  is  assumed  to  be  governed  by  the 
eigenvalues  of  the  Orr-Sommerf eld  equation  and  the  structure-borne  noise 
field  is  taken  to  be  flexural  in  nature.  The  structure-borne  noise  field 
is  decomposed  into  propagating  wave  fields  and  near-fields  such  as  would 
surround  foundation  structures  or  other  hull  impedance  discontinuities. 


It  is  found  that  for  the  frequency  range  generally  associated  with  struc- 
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ture-borne  noise,  say  30  - 10  Hz,  the  propagating  wave  fields  are  not 
destabilizing  in  that  they  do  not  contain  components  which  are  coinci- 
dent with  the  appropriate  eigenvalues  of  the  Orr-Sommerfeld  equation. 

By  contrast,  the  near-fields  will  contain  such  components.  These  near- 
field effects  are  interpreted  in  terms  of  equivalent  propagating  fields 
to  which  they  are  found  to  be  of  (relatively)  low  level.  Available  ex- 
perimental results  are  discussed  in  terms  of  the  submarine  environment 
and  it  is  found  that  although  transition  may  be  affected  by  structure- 
borne  noise  levels  less  than,  say,  the  expected  free-stream  turbulence 
levels,  the  measured  noise  levels  are  higher  than  one  expects  in  submarine 
hull  plating. 

Finally,  a mechanical  analogue  to  the  Orr-Sommerfeld  equation  is 
presented  in  the  form  of  a vibrating  elastic  plate  resting  on  a locally 
reacting  foundation  of  specific  form.  In  the  analogue  Reynolds  number 
becomes  the  cube  of  the  slenderness  ratio  of  the  plate.  Some  aspects 
of  the  solutions  to  the  Orr-Sommerfeld  equation,  as  it  is  applied  to 
boundary  layer  transition,  are  interpreted  in  terms  of  the  analogue 
system. 


I.  INTRODUCTION 

Early  studies  of  boundary  layer  control  were  concerned  primarily  with 

airfoils  where  suction  was  found  to  be  a promising  control  measure.^  Heating 

of  the  boundary  was  first  studied  for  gases,  for  which  a temperature  rise 

corresponds  to  an  increase  in  viscosity,  and  consequently  to  an  increase  in 

boundary  layer  thickness.  Heating  was  therefore  considered  as  a destabiliz- 
2 

ing  factor.  The  opposite  effect  is  achieved  in  water,  where  an  increase  in 
temperature  reduces  the  kinematic  viscosity  and  consequently,  the  boundary 
layer  thickness.  The  shift  of  emphasis  to  boundary  layer  control  in  water 
coincided  with  the  discovery  by  Kramer,  in  1960,  that  a compliant  boundary, 
specifically  a rubber  membrane,  can  delay  transition  from  laminar  to  tur- 
bulent flow.'* 

The  application  of  boundary  layer  control  techniques  to  submarines 
in  particular  is  desirable  since  boundary  layer  control  would  not  only 
reduce  drag,  but  also  turbulence-generated  self-noise  which  masks  sonar 
signals  at  the  higher  speeds.  Furthermore,  the  reduced  drag  entails 
reduced  propulsive  power  requirements.  This,  in  turn,  permits  a smaller 
hull  displacement,  reduced  auxiliary  power  requirements,  and  thus  a reduc- 
tion of  machinery-generated  radiated  noise  as  well.  Problems  that  may 
be  encountered  in  the  successful  application  of  boundary  layer  control 
techniques  to  submarines  concern  the  nature  of  the  environment  in  which 
the  submarine  operates.  That  is,  the  possible  destabilizing  effects  of 
free-stream  turbulence,  acoustic  disturbances,  hull  vibrations,  etc. 

Structure-borne  noise  has  no  effect  on  drag  in  current  submarines, 
because  most  of  the  hull  plating,  and  in  fact,  all  of  the  pressure  hull 
plating  not  concealed  by  thin  plating,  is  exposed  to  the  turbulent  boundary 
layer  at  all  speeds.  However,  if  the  boundary  layer  is  controlled  by,  say, 
heating  and  suction,  the  vibrating  pressure  hull  plating  will  be  in  contact 
with  the  laminar  boundary  layer.  This  report,  in  particular,  details  an 
investigation  into  the  potential  boundary  layer  destabilizing  action  of 
machinery  generated  structure-borne  noise  in  submarine  hull  plating, 
especially  the  effect  of  flexural  waves  in  submarine  hull  plating  on  tran- 
sition as  modelled  by  the  Orr-Sommerfeld  equation. 
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II.  the  characteristics  of  structure-borne  noise  in  submarine  hulls 


Machinery  generated  structure-borne  noise  in  operational  submarines 
consists  of  tonals  caused  by,  for  example,  blade  passage  in  pumps  as  well 
as  broad-band  components  caused  by,  for  example,  turbulent  flow  in  piping 
systems.  In  all  cases  the  excitation  levels  are  such  that  the  vibration 
response  of  the  hull  plating  is  the  result  of  a linear  phenomenon.  The 
spectrum  of  machinery  generated  structure-borne  noise  can  be  considered 
to  range  from  30  Hz  to  10  kHz. 

Machinery  induced  vibratory  forces  are  transmitted  through  a foundation 
structure  to  the  hull,  which  responds  to  these  excitations  primarily  in 
flexure.  Freely  propagating  flexural  waves  in  a plate  exhibit  dispersion. 
They  display  a phase  velocity,  in  vacuo 


c 
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J 

(c  hw)  2 
P 


(1) 
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where  c^  = compressional  wave  velocity  in  plates 
h = plate  thickness 
w = frequency  in  radian/sec. 

For  completeness,  it  should  be  mentioned  that  at  the  lower  frequencies 
this  result  is  modified  somewhat  by  the  effect  of  water  loading  which 
tends  to  reduce  the  velocity.  Also,  at  the  higher  frequencies,  the  above 
equation  must  be  modified  to  account  for  rotary  inertia  and  shear  effects. 
However,  both  corrections  are  not  significant  to  the  present  task. 

Hull  plating  typically  ranges  in  thickness  from  5 cm.  (2  in.)  for 
the  pressure  hull  to  1 cm.  (3/8  in.)  for  the  faired  portions  of  the  hull. 
The  overall  charcteristics  of  the  vibration  field  that  is  associated  with 
the  hull  plating  may  be  conveniently  summarized  by  plotting  the  appro- 
priate phase-velocity  vs.  frequency  relationship  for  relevant  parameters. 
This  is  shown  in  Fig.  1.  It  can  be  verified  that  the  flexural  wave 
velocity  in  the  pressure  hull  plating  exceeds  the  free  stream  flow  ve- 
locity over  the  entire  range  of  interest.  Furthermore,  it  is  noted  that 
the  flexural  wave  velocity  exceeds  the  sound  velocity  in  water  only  above 


the  coincidence  frequency,  which  is  approximately  5 kHz  for  representative 
5 cm.  (2  in.)  pressure  hull  plating. 


In  addition  to  freely  propagating  waves  in  the  hull  plating,  localized, 
or  near-field,  vibrations  are  present  in  the  immediate  vicinity  of  an  equip- 
ment-hull interface  as  well  as  at  any  major  discontinuity  in  the  impedance 
of  the  hull  structure.  At  low  frequencies  major  hull  discontinuities  are 
limited  to  bulkheads,  but  at  higher  frequencies  small  frames  also  provide 
significant  discontinuities  to  the  hull  plating.  The  most  significant 
difference  between  these  near-field  vibrations  and  propagating  waves  is 
that  for  a given  frequency  the  near-field  vibrations  contain  a continuous 
phase  velocity,  or  wavenumber,  spectrum.  If  one  models  this  effect  as  a 
uniform  line  excited  plate  of  infinite  extent,  or  equivalently  a point  ex- 
cited beam,  then  the  phase  velocity  spectrum  of  the  plating  response  at 


the  drive  point  is  as  shown  in  Fig.  2 for  a representative  value  of  struc- 
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tural  damping  given  by  n = 0.1.  Thus,  although  the  near-field  contains 


phase  velocity  components  c^,  such  that  c^/c^  >:>  1 > the  magnitudes  of  these 


components  are  low  relative  to  the  component  that  corresponds  to  the  char- 
acteristics of  freely  propagating  waves  in  the  plating,  i.e.,  c^/c  = 1. 


III.  BOUNDARY  LAYER  STABILITY;  TRANSITION 

The  classic  linear  analytical  approach  to  the  stability  of  two 
dimensional  parallel  flows  in  general,  and  boundary  layers  in  particular, 
is  to  solve  for  the  eigenvalues  of  the  Orr-Sommerf eld  equation.^  The 
Orr-Sommerfeld  equation  results  from  the  assumption  of  a disturbance  of 
the  form 


$ = <(>(y)  exp [ia (x  - ct)  ] 


where  x is  the  coordinate  along  the  direction  of  flow.  The  nature  of 
the  eigenvalues,  viz.,  (a,c),  is  then  used  to  determine  the  condition  of 
stability  of  the  flow,  i.e.,  whether  the  flow  can  support  a disturbance  of 
a given  form  without  producing  an  exponentially  amplified  response.  For 
example,  spatial  amplif ication  is  associated  with  negative  a^.  The  results 
of  these  analyses  are  typically  presented  in  the  form  of  stability  contours 
plotted  on  graphs  whose  abscissa  is  the  Reynolds  number  and  whose  ordinate 
is  the,  non-dimensionalized,  disturbance  wavenumber  for  various  phase  veloc- 
ities of  the  disturbance  (Fig.  3) . Alternatively,  contours  may  be  plotted 
on  graphs  of  frequency  w (~ac)  vs.  Reynolds  number  (Fig.  4).  The  neutral 
stability  contour  is  defined  as  the  contour  that  corresponds  to  a.  = 0. 
However,  observations  have  shown  that  in  many  instances  the  criterion 
of  neutral  stability,  i.e.,  that  transition  occurs  at  the  minimum  Reynolds 
number  capable  of  sustaining  a neutrally  stable  disturbance,  is  too  con- 
servative. Thus,  the  "en"  rule  has  been  developed. 

The  "e11"  rule  is  an  empirically  based  criterion  for  applying  linear 

6,  7 

stability  theory  to  the  problem  of  predicting  boundary  layer  transition. 

For  nearly  parallel  flow,  local  amplification  factors,  i.e.,  values  of  n, 
for  a given  disturbance  frequency  are  obtained  from  the  solutions  of  the 
Orr-Sommerfeld  equation  and  integrated  in  the  flow  direction  from  the  point 

g 

of  neutral  stability.  If  the  disturbance  whose  frequency  yields  the 
maximum  amplification  factor  at  transition  is  assumed  to,  in  fact,  cause 
transition,  then  it  is  found  that  the  value  of  the  amplification  factor 
itself,  i.e.,  the  quantity 
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Regardless  of  the  value  of  n,  the  point  is  that  transition  is  associ- 
ated with  disturbances  whose  wavenumber  and  phase-velocity  are  such  that 
they  coincide  with  the  neutral  stability  curves,  i.e.,  n = 0,  or  else  they 
penetrate  this  curve  to  some  degree,  i.e.,  n > 0.  This  fact  will  consti- 
tute the  basis  for  assessing  the  impact  of  hull  vibrations  on  transition. 


IV.  STRUCTURE-BORNE  NOISE  AS  A POTENTIAL  DESTABILIZING  AGENT 


Destabilization  of  a laminar  boundary  layer  by  structure-borne 
noise  requires  (1)  a transfer  mechanism  for  the  disturbance  to  interact 
with,  or  "penetrate,"  the  boundary  layer  and  (2)  once  penetrated,  an 
instability  mechanism  to  initiate  amplification  of  the  boundary  layer's 
response  to  the  disturbance. 

Four  different  configurations  through  which  structure-borne  noise 
in  submarine  hulls  represent  a potential  destabilizing  agent  are  shown 
in  Fig.  5.  Structure-borne  noise  in  uncoated  pressure  hull  plating  (Item 
1 in  Fig.  5)  represents  a potential  disturbance  which  is  introduced 
directly  into  the  boundary  layer.  The  drive  point  impedance  of  the  rela- 
tively thick  pressure  hull  plate  is  sufficiently  high  so  that  the  plating 
as  seen  by  the  boundary  layer  is  effectively  rigid.  In  other  words,  the 
pressure-hull  plating  can  be  assumed  to  exhibit  an  infinite  internal  im- 
pedance. Thus  the  relevant,  linearized,  analysis  requires  the  solution  to 
the  inhomogeneous  Orr-Sommerfeld  equation  with  the  inhomogeneity  being 
a prescribed  velocity  distribution  at  the  boundary.  The  importance  of 
this  circumstance  is  that,  although  the  actual  response  of  the  boundary 
layer,  e.g.,  amplification  factors,  etc.,  will  be  a function  of  the 
inhomogeneity,  the  instability  mechanisms,  i.e.,  the  eigenvalues  of  the 
Orr-Sommerfeld  equation,  correspond  to  the  homogeneous,  i.e.,  rigid 
boundary,  problem.  Thus,  the  criterion  for  destabilization  is  coinci- 
dence of  the  disturbance  wavenumber-phase  velocity  of  wavenumber- frequency 
characteristics  with  the  eigenvalues  corresponding  to  the  classical  Tollmien- 
Schlichting  instabilities  (see  Appendix  A) . 

A.  The  Effect  of  Propagating  Waves  on  Transition 

In  view  of  the  above  stated  requirement  for  destabilization  the 
question  as  to  whether  propagating  waves  in  submarine  plating  represent  a 
potential  destabilizing  factor  reduces  to  the  question  of  whether  these 
waves  are  coincident  with  the  eigenvalues  of  the  Orr-Sommerfeld  equation. 
Fortunately,  for  the  parameters  of  interest  coincidence  does  not  occur. 

This  can  be  seen  by  comparing  Figs.  1 and  3.  From  Fig.  3 it  can  be  seen 
that  Tollmien-Schlichting  instabilities,  for  a Blasius  profile,  exhibit 
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Q 

phase  velocities  which  are  bounded  by  c /U^  - 0.5.  In  order  to  compare 
the  results  of  Figs.  1 and  3,  it  is  convenient  to  express  this  velocity 
as  follows: 

CTS  " 0,5  U°°  = °-5caM  = °-5  cf(ca/cf)M  = [0.5  M (^/to)  2 ] cf  (2) 

1 2 

where  M is  the  Mach  number  and  u (=  12  c /(he  ))  is  the  critical 

c p 

frequency  of  the  hull  plating  in  water.  Thus,  for  coincidence  between 
the  phase  velocity  of  the  flexural  waves  in  the  plating  with  that  of  T-S 
waves 


to  - (0.5M)  ^co 

c 

For  a speed  of  80  knots,  this  requirement  for  the  pressure  hull  plating 
implies  a frequency  less  than  1 Hz,  which  is  outside  our  frequency  range 
of  interest.  For  completeness  it  is  noted  that  this  conclusion  is  not 
changed  if,  instead  of  the  Blasius  profile,  one  uses  the  profiles  expected 
for  heated  boundary  layers  as  calculated  by  Rockwell.10 

Item  3 in  Fig.  5 poses  a different  circumstance  in  that  here  the 

structure-borne  noise  field  does  not  interface  with  the  boundary  layer, 

but  rather  is  isolated  by  means  of  a compliant,  anti-radiation,  coating. 

This  differs  from  the  previous  situation  since  here  the  boundary/boundary- 

layer  interface,  by  design,  exhibits  a relatively  low  impedance.  This 

means  that  a feedback  loop  is  indeed  possible  between  the  boundary  layer 

and  the  boundary.  In  other  words  the  non-rigid  boundary  can  modify  the 

eigenvalues  of  the  Orr-Sommerfeld  e~  on  and  in  fact  introduces  additional 
9 

sets  of  eigenvalues.  It  must  therefore  be  determined  if  these  modified 
eigenvalues,  i.e.,  instabilities,  can  be  excited  by  the  structure-borne 
noise  field  associated  with  the  pressure  hull  plating.  An  analysis  of 
this  situation  is  presented  in  Appendix  B where  it  is  again  shown  that  for 
the  relevant  parameters  the  conditions  for  destabilization  are  not  satisfied 
by  present  day  submarine  constructions. 

Item  2 ir.  Fig.  5 falls  somewhere  between  Items  1 and  3.  However, 
in  lieu  of  the  previous  results  it  is  concluded  that  structure-borne  noise 
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in  the  faired  hull  plating,  in  the  form  of  propagating  waves,  as  with  the 
pressure  hull  plating,  is  not  a potential  destabilizing  agent.  Finally, 
for  completeness,  disturbances  associated  with  Item  4 in  Fig.  5,  propa- 


gating acoustic  waves  of  moderate  level  also  fall  outside  the  wavenumber/ 
phase-velocity  range  of  potential  instabilities  as  calculated  via  the  Orr- 
Sommerfeld  equation. 


It  should  be  noted  that  in  addition  to  the  mechanisms  of  modifying 
boundary  layer  stability  indicated  in  Fig.  5,  a mechanism  based  on 


vortex  generation  by  an  acoustic  field  is  suggested  by  the  work  of  Mechel 
and  Schilz.^  This  mechanism  requires  a boundary  layer  flow  along  a sur- 


face that  has  discontinuities  in  the  form  of  sharp  edges.  The  oscillating 
flow  associated  with  an  intense  sound  field  can  establish  an  alternating 


vortex  structure  when  the  acoustic  motion  is  normal  to  a sharp  edge.  The 


vortex  motion  then  couples  in  a nonlinear  manner  with  the  mean  flow  to 


cause  higher  frequency  components  of  the  disturbance  to  be  generated. 

The  possibility  therefore  can  not  be  dismissed  that  the  intense  acoustic 


motion  results  in  flow  disturbances  that  can  match  either  constructively 


or  destructively  the  incipient  Tollmien-Schlichting  disturbances.  It 


should  be  noted  that  this  sharp  edge  mechanism  would  not  be  expected  to 


be  present  on  a smoothly  faired  shape  such  as  that  of  a submarine. 


B.  The  Consequences  of  Near-Field  Effects  on  Transition 


As  discussed  earlier,  near-field  vibrations  differ  from  the 


propagating  field  in  that  the  near-field  contains  a continuous  wave- 


number, or  phase  velocity,  spectrum.  This  means  that  the  near-field 
will,  in  fact,  contain  a component  that  is  coincident  with  T-S  waves. 
Therefore,  here  the  question  as  to  whether  near-field  effects  in  sub- 
marine plating  represent  a potential  destabilizing  factor  reduces  to  a 
question  of  level.  For  hull  plating  in  the  immediate  vicinity  of  an 
equipment  foundation  or  for  substantial  impedance  discontinuities,  the 


overall  near-field,  i.e.,  non-propagating  levels,  can  be  expected  to  be 


of  the  same  level  as  the  propagating  wave  levels.  The  relative  magnitudes 
of  those  phase  velocity  components  of  the  near-field  that  may  be  coinci- 


dent with  T-S  waves  can  be  obtained  from  Fig.  2.  Note  that  the  excitation 


spectrum  level  decreases  monotonically  with  a decrease  in  phase  velocity 
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[ 

1 


from  the  point  corresponding  to  the  phase  velocity  of  propagating  waves. 

A measure  of  the  potential  destabilizing  effect  of  this  near-field 
may  be  obtained  by  relating  it  to  an  "equivalent"  propagating  field. 

The  meaning  of  equivalent  is  equal  wavenumber  and  phase  velocity.  As  an 
example  consider  the  flexural  response  of  thin  hull  plating  at  the  high 
end  of  the  frequency  range  of  interest,  say  10  kHz.  From  Eq.  2,  at  this 
frequency  and  assuming  a speed  of  80  knots,  coincidence  occurs  at 


c /c,  = c„c/cf  : 2 x 10  From  this  value  the  spectrum  level  continues 
Y t lb  r 4 

to  decay  at  a rate  proportional  to  (c  /cf)  . Referring  to  Fig.  2 the 

—8 

level  of  this  component  is  down  by  a factor  of  1.6  x 10  from  the  level 


of  the  propagating  component.  Thus,  the  near-fields  surrounding  local 
excitations  or  discontinuities  represent  a comparable  threat  to  stabiliza- 
tion as  the  equivalent  propagating  field  whose  level  is  down  by  a factor 
of  1.6  x 10-8  or  156  dB. 


r 


V.  EXPERIMENTAL  DATA  RELEVANT  TO  THE  PROBLEM  OF  STRUCTURE-BORNE  NOISE 


EFFECTS  ON  TRANSITION 

} 

A.  Destabilization  with  Coincidence 

In  this  section  the  salient  features  of  an  experimental  study 

of  the  effects  of  boundary  vibrations  on  transition  are  presented.  The 

12 

work  was  performed  by  Schilz  in  1969.  Even  though  this  work  was 
funded  by  the  U.S.A.F.  European  Office  of  Aerospace  Research  and  published, 
albeit  in  German,  it  appears  to  have  attracted  little  notice  in  the  U.S. 

1 . Description  of  Tests 

It  will  be  found  in  these  tests  that,  as  discussed  in 
Section  IV.,  the  effects  of  flexural  waves  propagating  in  a plane  boundary 
are  most  significant  if  the  disturbance  is  coincident  with  a T-S  instability. 

The  T-S  phase  velocity  lies  between  1 and  20  m/s  in  Schilz's  tests.  No 
practical  structural  material  can  achieve  these  low  velocities  in  the  fre- 
quency range  of  interest.  Schilz  overcame  this  difficulty  by  using  a 

j 

transducer  array  wherein  adjoining  elements  are  driven  out-of-phase.  One 
pair  of  elements  thus  simulates  one  wavelength.  Since  the  spacing  of  the 
transducers  in  a given  test  is  invariant,  and  can  only  be  varied  discon- 
tinuously  between  tests,  increasing  values  of  Um  are  associated  in  each 
test  with  rising  frequencies.  The  boundary  condition  thus  achieved  dis- 
plays, not  surprisingly,  a strong  harmonic  content. 

The  displacement  amplitudes  of  the  vibrating  boundary  were  measured 
with  a capacitance  probe.  Hot  wire  anemometers  were  used  to  measure  the 
flow  velocities.  The  tests  were  performed  on  a flat  plate  with  a sharp 
leading  edge.  Unless  tripped,  the  laminar  boundary  layer  extends  over  the 
entire  plate  (0.5  m) . 

2 . Experimental  Results 

Destabilization  occurs  primarily  within  the  neutral  sta- 
bility curve  of  T-S  waves,  but  was  occasionally  observed  at  lower  fre- 
quencies. Only  semiquantitative  measurements  could  be  obtained  for  the 
threshold  amplitude  of  flexural  waves  required  to  trip  the  boundary 
layer  for  different  values  of  U^.  It  was  found  that  for  the  coincidence 
conditions,  the  amplitudes  required  to  affect  the  boundary  layer  are 
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much  smaller  (10  millimeters  or  less)  than  the  boundary  layer  displace- 
ment thickness  (0.4  millimeters).  The  results  show  a marked  decrease  with 
increasing  U^. 

A marked  rise  of  the  threshold  occurs  on  both  sides  of  the  coincidence 
frequency,  defined  as  the  frequency  for  which  the  effective  phase  velocity 
of  the  "flexural  wave"  matches  the  T-S  velocity.  When  the  disturbance  is 
in  the  form  of  a band  of  noise,  only  that  band  which  contains  the  coin- 
cidence frequency  shifts  the  transition  point  upstream. 

3.  Conclusions 

The  principal  outcome  of  the  experiments  is  that  the  con- 
clusion drawn  in  Section  IV.,  namely  that  it  is  improbable  that  machinery 
generated  structure-borne  noise  in  the  form  of  propagating  waves  could 
destabilize  a laminar  boundary  layer,  appears  to  be  verified. 

Another  conclusion,  drawn  by  Schilz  and  not  unrelated  to  the  problem 
at  hand,  is  that  structure-borne  noise  may  provide  a means  for  actually 
suppressing  turbulence,  i.e.,  delaying  transition.  For  this  hypothesis, 
Schilz  relies  heavily  on  Benjamin's  analysis.^  Benjamin  defines  a type 
A instability  as  the  T-S  wave  modified  by  the  presence  of  a flexible 
boundary.  He  shows  that  this  is  an  energy-deficient  wave  mode.  Damping 
at  a flexible  boundary  further  enhances  this  energy  deficiency,  thus 
leading  to  energy  flow  from  the  main  stream  to  the  boundary  layer,  and 
consequently  further  destabilization.  Schilz  postulates  that  instead  of 
supplying  the  needed  energy  by  diverting  it  from  the  main  stream,  as  sug- 
gested by  Benjamin,  this  energy  can  be  provided  by  actively  exciting 
flexural  waves  in  the  boundary.  These  waves  convey  energy  to  the  boun- 
dary layer,  thus  compensating  for  energy  dissipation  in  the  T-S  waves. 

The  flexural  waves  thus  play  the  role  of  negative  damping,  which  can  be 
anticipated  to  have  a stabilizing  effect,  since  positive  damping  exerts 

a destabilizing  influence.  At  coincidence,  the  condition  for  energy 

14 

flow  from  the  wall  into  the  fluid  depends  on  the  phase.  (Becker  (includes 
that  for  any  stable  inviscid  velocity  profile,  energy  flow  is  from  the 
fluid  to  the  vibrating  boundary.  This  can  presumably  be  assimilated  to 
positive  damping  by  the  boundary,  which  Benjamin  showed  to  be  destabilizing. 
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It  would  appear  the  Schilz's  technique  of  forcing  energy  to  flow  from  the 
wall  into  the  boundary  layer  can  only  be  implemented  for  unstable  profiles.) 

Schilz  finds  that  he  does  indeed  achieve  stabilization  over  a broad 
phase  angle  A4>  centered  approximately  on  180°,  i.e.,  for  phase  opposition 
between  T-S  waves  and  the  wall  vibrations.  These  results  were  obtained 
as  follows:  the  boundary  layer  is  tripped  to  generate  T-S  waves,  the  re- 
sulting transition  point  being  located  downstream  of  the  flexural  array. 

The  array  is  then  energized  to  produce  flexural  waves  whose  velocity 
coincides  with  that  of  the  T-S  waves.  Appropriate  amplitudes  and  A<J>  : 180° 
can  lead  to  stabilization.  For  a 180°  phase  angle,  an  increase  in  the 

flexural  wave  amplitude  first  increases  the  critical  Reynolds  number  and, 

_ ^ 

for  amplitudes  exceeding  an  optimal  amplitude  3 x 10  meters  for  = 

27.4  m/s,  reduces  it.  If  the  boundary  layer  is  not  tripped,  these  extremely 
low  flexural  wave  amplitudes  have  a destabilizing  effect,  presumably  be- 
cause they  are  in  phase  with  the  T-S  waves  they  generate,  rather  than  being 
out  of  phase  with  the  T-S  waves  generated  by  the  tripping  wire. 

B.  Destabilization  Without  Coincidence 

In  a separate  paper  Schilz  discusses  non-linear  effects  that 
may  be  important  when  coincidence  is  not  achieved.^0  A more  detailed 
description  of  this  type  of  circumstance,  i.e.,  non-coincidence,  can  be 
found  in  the  reports  by  Norair,^  documenting  their  preliminary  wind 
tunnel  experiments  on  the  effects  of  noise  on  transition. 

1 . Description  of  Tests 

A 4-percent-thick  straight  laminar  suction  wing  of  17- 
foot  chord  was  investigated  in  the  Norair  7-  by  10-foot  low  turbulence 
wind  tunnel  at  a = 0°  angle  of  attack  in  the  presence  of  external  sound, 


and  in  addition,  panel  vibration.  The  sound  consisted  of 


discrete 


frequencies  and  octave  bands  of  random  noise  in  the  150  to  4000  cps 
frequency  range,  while  the  vibration  frequencies  were  100,  190,  and 
1290  cps. 

Considerations  of  the  type  of  structure  contemplated  by  Norair 
indicated  that  peak  panel  accelerations  of  about  lOg  in  the  100  to  1000 
cycle  frequency  range  could  be  expected.  Thus,  to  answer  the  question  as 
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to  the  potential  influence  of  mechanical  vibrations  of  an  external  wing 
surface  on  the  behavior  of  a laminar  (suction)  wing,  a high-speed,  d.c. 
motor-driven,  reciprocating  mass  shaker  was  attached  to  the  wing.  The 
motor  power- limited  the  shaker  to  a frequency  of  about  260  Hz  and  a 46  g 
peak  panel  acceleration  for  intermittent  operation. 

2.  Experimental  Results 

In  all  instances  of  lOg  peak  acceleration  or  more,  suction 
had  to  be  increased  above  the  minimum  drag  value  to  maintain  full  chord 
laminar  flow.  For  a peak  panel  velocity  ratio  v/Uk  = 2.0  x 10  \ the 
maintenance  of  full  chord  laminar  flow  by  raising  suction  level  required 
a 24  percent  increase  in  suction  coefficient  and  a corresponding  5 percent 
increase  in  equivalent  drag,  compared  to  a case  with  no  vibration.  Since 
panel  velocity  ratios,  v/U I , in  the  order  of  2 x 10  3 are  possible  for 
thin  straight  wing  laminar  suction  aircraft,  it  was  concluded  that  some 
additional  suction  might  be  required. 

To  interpret  the  magnitude  of  the  vibration  levels  in  the  Norair 
tests  to  the  problem  of  structure-borne  noise  in  submarine  hull  plating, 
consider  the  example  of  a machine  with,  say,  a motor  imbalance  that  pro- 
duces a 1 lb.  vibratory  force  on  the  hull  plating.  If  one  crudely  models 
the  hull  plating  as  an  infinite  elastic  plate  in  order  to  obtain  an  es- 
timate of  its  impedance,  then  2 inch  plating  yields  a drive  point  imped- 

3 

ance  z - 10  lb. /in. /sec.  Thus,  the  excitation  will  produce  hull  accelera- 
tions of  approximately 

a/g  = wv/g  = 2n  (1)  (10~3/386) f 
- 1.5  x 10~5  f 

where  f is  the  excitation  frequency  in  Hz.  If  one  considers  5 kHz  the 
highest  frequency  of  interest  the  maximum  plating  acceleration  levels 
are  given  by  a/g  * 10  3 compared  to  a/g  » 10  in  the  Norair  tests. 

This  particular  result  is  somewhat  disturbing  in  that  although 
destabilization  required  structure-borne  noise  levels  two  orders  of 
magnitude  greater  than  one  expects  in  submarine  hull  plating,  the 
mechanism  behind  the  destabilization  process  appears  to  be  relevant. 
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For  example,  the  structure-borne  noise  levels  in  the  Norair  test,  as  well 
as  in  an  operational  submarine,  are  small  in  terms  of  both  the  boundary 
layer  thickness  and  the  characteristic  structural  dimensions,  viz.  plating 
thickness  and  length.  Unforturnately,  at  present  the  available  Norair 
documentation  covers  their  preliminary  results.  If,  in  fact,  their  work 
or  similar  work  is  continuing,  it  would  be  most  relevant  to  the  problem 
at  hand. 


VI.  A MECHANICAL  ANALOGUE  TO  BOUNDARY  LAYER  STABILITY 


A.  The  Analogue  System 

During  the  course  of  this  study  a mechanical  analogue  to  the 
Orr-Sommerf eld  equation  was  discovered.  This  analogue  was  pursued  and 
the  results  are  presented  in  this  section.  A qualitative  description  of 
some  features  of  the  submarine  structure-borne  noise  problem  is  presented, 
i.e.,  interpreted,  in  terms  of  the  analogue  system  in  Section  B. 

Consider  the  steady  state  vibrations  of  an  elastic  plate.  For  present 
purposes,  the  plate  is  taken  to  be  of  infinite  extent  along  the  x direction 
of  "characteristic"  width  f and  thickness  h (Fig.  6) . The  partial  differ- 
ential equation  governing  the  vibration  response  of  a plate  can  be  obtained 
from  the  Timoshenko-Mindlin  model  of  plate  vibrations.'*’  This  model  allows 
for  the  same  structural  wavelength,  two  different  modes  of  wave  propagation 
corresponding  to  two  different  ratios  of  shear  to  bending  motion.  Or, 
stated  differently,  for  the  same  frequency  this  model  allows  two  different 
modes  of  wave  propagation  corresponding  to  two  different  structural  wave- 
lengths. The  model  is  generally  considered  valid  up  to  frequencies  that 

produce  structural  wavelengths  A such  that  (A  /h)  - 0(10). 

s s 

Using  the  rotation  of  the  plating  cross-section  in  the  y-z  plane, 
Y(x,y),  as  the  response  parameter  of  the  plate,  the  governing  partial 
differential  equation  of  motion  may  be  expressed  as 

[V2  + a2][V2  + (cd/cs)  V]y(x,y)  - 82{(22Y(x,y) 

+ 8 IV2  + (c,/c  ) 2fi2)M (x,y) } = 0 
d s 

where 

2 2 2 2 2 

= ravaxz  + a /ay  ) 

x,y  = coordinate  axes  nondimensionalized  to  H 

l = characteristic  plate  dimension  along  the  y axis 
h = plate  thickness 

6 = slenderness  ratio  of  plate  = / 12 ( l/h) 
c^/cs  = ratio  of  dilatational  to  shear  wave  speed  in  plate 


(3) 
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ft  = nondimensionalized  frequency  = uii/c^  where  u>  is  the  circular 
frequency 

Q(x,y)  - distributed  lateral  force  acting  on  the  plate,  nondimensionalized 
to  12E/U  - v2) 

M(x,y)  = distributed  moment  acting  on  the  plate,  nondimensionalized  to 
12KH/(1  - v2) 

b , v = elastic  modulus  and  Poisson's  ratio  of  plating  material 


Now  consider  the  case  where  the  lateral  forces  and  moments,  Q and  M, 
represent  the  response  of  a homogeneous  and  isotropic,  but  locally 
reacting,  foundation  on  which  the  plate  is  resting  (Fig.  6) . This  may 
be  expressed  in  terms  of  the  foundation's  impedances  Z^(y)  and  ZM(y) 


Q(x,y)  = -iftz^(y) y(x,y) 

M(x,y)  = -iftZ  (y)y(x,y) 
M 


(4) 


Further,  the  assumption  is  made  that  for  the  plate  under  consideration 

<c,/c  )2  = 1 (5) 

d s 

Althouqh  this  is  not  ordinarily  the  case  for  homogeneous  plating  materials, 
one  can  envision  this  being  accomplished  by  means  of  a sandwich-like  plate 
construction. 

Substituting  Eqs.  4 and  5 into  Eq.  3 yields 

(V2  + n2)(V2  + ft2 J y (x, y)  - e2{fi2r(x,y)  - ift0fafZQ(y)Y(x,y))/3y 

+ (V2  + ft2)  iZ  (y)  y (x,y)  ] ] } = 0 (6) 

M 

It  follows  that  the  wave  number  spectrum  in  the  x direction  oi  the  plate's 
response,  as  given  by  the  Fourier  transform  of  Eq.  6,  satisfies  the 
equation 
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2 " 4 2 2 1 

y (y)  - 2n  y (y)  + n y<y)  = 3 (ft  Y(y)  - ii2B  [ ( zq < y) y <y)  I 

+ [2M(y)Y(y)l  - n2[ZM(y)Y(y)]"U 


(7) 


where 


y(y) 


Y (x,y) expdk^s)  dx. 


2 M 2 n2, 
n = (k  - 52  ) 
x 


and  the  "prime"  notation  is  used  to  represent  differentiation  with  respect 
to  y.  At  this  point  the  form  of  the  foundation  impedances  are  restricted 
to  take  the  following  form: 


(i)  Z (y)  = -2ZM(y) 


(8) 


and 


(ii)  e|ZM(y)  + n ZM(y)  I >>  52 


Thus  Eq.  7 becomes 


"«  2 " 4 3 2 

Y (y)  - 2n  y (y)  + n Y<y)  = 152B  Un  z (y) 


+ zu(y) lY(y)  - zM(y)Y  (y) 

M M 


or 


""  2 " 4 3,2*" 

y (y)  - 2n  y (y)  + n y (y)  = -inB  (In  ZM  (y)]Y(y) 


- zM(y)y  (y) 


where 


Z*(Y)  = -zM(y)/t(kx/52)2  - 11 i 
= -(52/n)zM(y) 


(9) 


(10) 


(11) 


(12) 
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The  solutions  to  Eq.  11  yield  the  eigenvalues  (n)  as  a function  of  £5 


and  ZM(y) . 

Eq.  11  will  now  be  compared  to  the  Orr-Sommerfeld  equation  which, 
for  completeness,  is  written  below 

2 " 4 2 

♦ (y)  — 2at  <J>  (y)  + a 4>  (y)  = -iaR{[(U(y)  - c)a 

+ U (y)  ] <t>  (y)  - (U(y)  - c ) <f>  (y) } (13) 


<}> (y)  is  defined  by  the  perturbation  stream  function  4>(x,y)  which  is  assumed 
to  be  of  the  form 

4,(x,y)  = <(>  (y)  exp]  ia  (x  - ct)  (14) 

x and  y are  cartesian  coordinates  parallel  and  normal  to  the  direction  of 
flow,  U(y)  describes  the  main  flow,  and  R is  the  Reynolds  number. 

In  dealing  with  boundary  layers,  U(y)  represents  the  mean  velocity 
profile  within  the  boundary  layer  and  all  equations  can  be  non-dimen- 
sionalized  whereby  lengths  are  normalized  to  the  boundary  layer  thickness 

5 ; velocities  to  the  free  stream  velocity  U , time  by  5/U^  and  pressure 

2 °° 

by  pU^  where  p is  the  fluid  density. 

For  a given  main  flow  profile  the  eigenvalues  of  the  Orr-Sommerfeld 
equation,  as  a function  of  Reynolds  number,  are  given  by  the  wavenumber 
and  phase  velocity  set  (a,c)  that  satisfies  Eq.  13.  Referring  to  Eq.  14 
it  is  noted  that  < 0 denotes  spatially  amplified  disturbances  with 
increasing  values  of  x.  The  subscript  "i"  is  used  to  denote  the  imaginary 
component. 


A comparison  between  Eqs.  12  and  13  reveals  the  following  formal 
analogy  (Table  1) . 


Table  1 - Mechanical  Analogue  to  Parallel 
Flow  Stability 


Parallel  Flow  Stability 
(Orr-Sommerf eld  Model) 


Plate  Vibrations 
(Timoshenko-Mindlin  Model) 


It  is  interesting  to  note  the  following: 

1.  Reynolds  number  (R)  is  analogous  to  the  cube  of  the  slenderness 
ratio  of  the  plate  (3)  whereby  large  Reynolds  number  corresponds  to  slender 
plating. 

2.  The  inviscid  components  of  the  Orr-Sommerf eld  equation  are 
modelled  by  the  plating  while  the  mean-flow  related,  or  viscid,  components 
are  modelled  by  the  locally  reacting  foundation  upon  which  the  plating 
rests . 


3.  The  critical  layer,  that  is,  the  value  of  y across  which  the 
real  component  of  the  quantity  [U(y)  - c]  changes  sign,  corresponds  to 

the  value  of  y in  the  analogue  across  which  the  real  component  of  the 

* 

foundation  impedance  lzM(y) 1 changes  sign. 

4.  Neutral  stability,  whereby  a and  c are  purely  real  quantities, 
corresponds  to  a completely  conservative  analogue  system  with  the  stipu- 
lation that  k >12. 

x 
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To  complete  the  formal  analogy  the  appropriate  boundary  conditions 
must  be  included.  This  appears  to  be  a potential  strength  of  the  analogy 
since  a vibrating  plate  can,  in  general,  admit  passive,  as  well  as  active, 
and  conservative  as  well  as  dissipative,  terminations.  Thus,  for  example 
parallel  flows  between  general  impedance  boundaries  should  be  amenable 
to  the  analogue.  However,  for  present  purposes,  it  is  merely  noted  that 
the  classical  boundary  conditions  for  boundary  layer  flow  over  an  effec- 
tively rigid  plane,  which  are  given  by 

4>  (0)  = 4>  ’ (0)  =0  (15a) 


(15b) 

are  modelled  in  the  analogue  by  a plate  that  is  semi-infinite  in  extent 
along  the  y axis  and  cantilevered  from  an  effectively  rigid  termination 
at  the  y = 0 boundary,  i.e., 

y(k  ,0)  = y’ (k  ,0)  = 0 (16a) 

' x x 

and 

I 

(16b) 

Thus  in  theory,  although  it  is  not  necessarily  proposed  here,  if 
one  were  to  construct  such  a structure  as  is  described  above,  and  impulse 
excite  it  such  that  all  frequencies  and  wavenumbers  were  excited,  then 
the  plate's  response  would  contain  the  solutions  to,  i.e.,  the  eigen- 
values of,  the  Orr-Sommerfeld  equation. 


lim  y (k  , y)  = lim  Y(k  ,y)  = 0 

X 

y-x»  y-*oo 


and 


lim  4>  (y) 

y-x» 


lim  <j>  (y)  = 0 

y-x» 
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B.  Stability  Phenomena  in  the  Analogue  System 


In  the  previous  section  a formal  analogue  to  the  Orr-Sommerfeld 
equation  was  presented.  In  this  section  some  features  of  the  solutions 
to  the  Orr-Sommerfeld  equation,  when  applied  to  boundary  layers,  are  in- 
terpreted in  terms  of  the  analogue  system.  However,  it  is  noted  that 
any  application  of  these  interpretations  to  experimental  observations  of 
transition  presupposes  that  the  mechanisms  governing  the  observations 
are,  in  fact,  captured  by  the  linearized  Orr-Sommerfeld  theory.  It  is 
also  noted  that  in  this  section  the  spatial  amplification  formulation 
to  the  solutions  of  the  Orr-Sommerfeld  equation  will  be  used.  That  is, 
the  wavenumber  will  be  considered  to  be  complex,  a = + ia^,  and  fre- 
quency will  be  taken  to  be  real  ui  = ac  = This  leads  to  a spatial 

amplif ication  factor  of  the  form  exp[-ot^x)  and  thus  negative  values  of 
yield  amplification  and  positive  values  of  correspond  to  attenu- 
ation. 

1.  Reduced  Analogue  for  Divergence 

"Divergence"  as  opposed  to  "flutter"  is  the  particular 
form  of  instability  whereby  wave  amplitude  varies  monotonically  with 
position.  Divergence  corresponds  to  solutions  of  the  Orr-Sommerfeld 
equation  with  purely  imaginary  values  of  a,  i.e.,  a = ia^.  In  the 
analogue  system  this  implies  that  n is  purely  imaginary  and  thus 
(k^/ft)  <1.  In  fact  it  is  of  particular  interest  to  consider  that  this 
condition  be  satisfied  in  the  limit  (k^/ft)  = 0.  Limiting  ourselves  to 
this  situation  the  corresponding  reduction  in  the  analogue  is  shown  in 
Table  2. 

To  interpret  this  reduced  analogue  consider  the  analogue  system  in 
Table  1 as  a two-dimensional  structural  waveguide,  albeit  leaky  along  the 

y axis.  Since  the  plating  itself  is  uniform  the  extent  of  this  "leak" 

* 

will  depend  on  the  variations  in  ZM(y),  with  strong  variations  tending 
to  reflect  rather  than  transmit  energy  and  thus  minimize  the  leakage. 

Each  vibration  mode  of  such  a system,  i.e.,  a mode  with  a given  mode  shape 
in  the  y direction,  will  generally  exhibit  a high-pass  cut-off  frequency, 
i.e.,  a frequency  below  which  the  mode  cannot  freely  propagate  along  the 
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Table  2 - Reduced  Analogue  for  Boundary 
Layer  Divergence 

— 

Boundary  Layer  Stability  Beam  Vibr. 


1 

Boundary  Layer  Stability 

Beam  Vibrations 

(Orr-Sommerfeld  Model) 

(Timoshenko-Mindlin  Model) 

6 

l 

' 

y 

y 

<t>  (y) 

Y (0,y) 

- 

+ ifi 

e3 

| U(y)-c 

±Z*(y)=±iZM(y) 

x-direction.  This  frequency  will  be  the  natural  frequency  of  the  corre- 

f . , TO 

sponding  standing  wave  mode  of  the  waveguide  cross-section  along  the  y axis, 

| and  thus  may  be  calculated  by  setting  the  x-direction  wavenumber  (k  ) equal 

to  zero  in  the  equation  of  motion  (Eq.  11) . Therefore,  the  wavenumber  (a) 

of  the  divergent  instability  associated  with  a given  phase  velocity  (c) 

may  be  interpreted  as  the  cut-off  frequency  (Q)  of  the  analogue  system 

associated  with  a given  foundation  impedance  Z (y) . The  reduced  analogue 

M 

shown  in  Table  2 contains  a number  of  interesting  features: 

1.  Physically,  the  dimensionality  of  the  analogue  has  been  reduced 

I to  that  of  a beam  (in  the  y-z  plane)  rather  than  a plate. 

2.  Real  values  of  c imply  a purely  imaginary  foundation  impedance 
ZM(y) . This  implies  a non-dissipative  foundation  and  thus  the  mechanical 
analogue  is  a purely  conservative  system. 

3.  Assuming  frequency  to  be  a real  positive  quantity,  the  positive 
M 

sign  for  Z (y)  implies  a = +ifi,  and  thus  in  view  of  Eq.  14,  corresponds 

M 

to  stability.  Similarly,  the  negative  option  for  Z (y)  corresponds  to 
an  analogy  where  a = -ifi  and  yields  the  unstable,  viz.  divergent,  regime. 


. Jl 
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2.  Transition  in  the  Reduced  Analogue  System 

As  discussed  in  Section  III,  the  e°  rule  is  an  empirically 
based  rule  whereby  if  the  disturbance  whose  frequency  yields  the  maximum 
amplification  factor  at  transition  is  assumed  to  cause  transitioi  , then 
it  is  found  that  the  value  of  the  amplification  factor  itself  is  reasonably 
constant  and  corresponds  to  n - 9.  This  result  is  extraordinary  in  that 
it  is  difficult  to  justify  the  linear  theory  in  connection  with  such  large 
amplification  rates.  Now  consider  this  criterion  in  the  reduced  analogy 
system  for  divergence  in  the  case  of  truly  parallel  flow  where  the  am- 
plification factor  is  constant  in  the  flow  direction  and  given  by 

a.dx  = -a. Ax  = n 
1 l 

where  the  subscripts  N.S.  and  T refer  to  neutral  stability  and  transition  respectively. 

From  Table  2 an  instability  corresponding  to  a given  amplification  factor 

is  analogous  to  the  response  of  the  mechanical  system  at  the  frequency  given 

by  (2  = ia  = i(ia^)  = - a ^ = n/Ax.  In  other  words,  in  the  analogue  system 

the  en  rule  is  not  to  be  interpreted  in  terms  of  a critical  amplification 

factor,  but  rather  in  terms  of  a critical  frequency. 

3.  Coupling  of  Disturbances  to  the  Boundary  Layer 
in  the  Reduced  Analogue  System 

Generally  transition  is  assumed  to  require  both  wavenumber 
and  frequency  (or  phase  velocity)  matching  in  the  flow  direction  with  a 
potential  disturbance.  However,  since  the  analogue  in  Table  2 has  been 
reduced  to  the  zero  x-direction  wavenumber  component  of  the  response  of 
the  plating  structure  by  setting  wavenumber  matching  between  the 

structure  and  a potential  disturbance  merely  requires  that  the  disturbance 
exhibit  a zero  x-direction  wavenumber  component  in  its  spectrum.  In  fact, 
since  the  dimensionality  of  the  reduced  analogue  is  that  of  a beam  (in 
the  y-z  plane)  rather  than  a plate,  the  condition  of  stability  is  a purely 
local  phenomenon  with  no  coupling  in  the  x,  or  flow,  direction.  A result- 
ing instability  might  well  be  called  a "turbulent  spot." 
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The  eigenvalues  of  the  Orr-Sommerfeld  equation  that  are  modelled 
by  this  reduced  analogue  may  correspond  to  the  category  of  disturbance 

generally  referred  to  as  Class  C,  since  for  this  class  there  is  apparently 

9 

no  dependence  of  the  amplification  rates  on  wavenumber. 

It  is  also  of  interest  to  consider  the  behavior  of  the  reduced 
analogue  structure  for  large  values  of  y.  Here  the  foundation  impedance 
becomes  uniform  since  U(y)  approaches  the  free  stream  velocity  U^.  The 
resulting  structure,  that  of  a uniform  beam  resting  on  a uniform  foundation, 
will  exhibit  a cut-off  along  the  y axis,  i.e.,  normal  to  the  boundary  layer, 
if  the  foundation  is  stiffness-like.  This  frequency  may  be  calculated  from 
Eq.  11  by  setting  k^,  the  y-direction  wavenumber  component  of  the  solution, 
equal  to  zero  as  well  as  k^  = 0.  This  yields  the  cut-off  frequency  for 
large  y 

^cut-off  = “6  ZM 

★ 

with  the  requirement  that  ZM  be  negative.  From  Table  2 it  is  seen  that 
the  analogue  for  divergence  requires  that  ^ = -d  - c)  and  thus 

for  potential  disturbances  with  phase  velocities  typically  associated 

* 

with  transition,  viz.  c < 1,  Z„  is  indeed  negative.  Therefore  for  a 

M 

divergent  instability  with  c < 1 a cut-off  frequency  exists  and  it 
increases  with  increasing  slenderness  ratio,  i.e.,  Reynolds  number,  such 
that  for  practical  parameters  this  form  of  instability  can  only  be  excited 
by  a disturbance  in  the  near-field  of  the  boundary  layer. 

4 . Submarine  Structure-Borne  Noise  in  the  Analogue  System 

Throughout  this  report,  and  indeed  throughout  most  of  the 
relevant  published  literature,  destabilization  is  assumed  to  require  a 
coincidence  between  the  disturbance  and  those  eigenvalues  of  the  Orr- 
Sommerfeld  equation  that  imply  wave  amplif ication.  In  other  words,  for 
example,  the  possibility  of  a phase  velocity  or  wavenumber  conversion 
mechanism  within,  or  in  the  immediate  vicinity  of,  the  boundary  layer  is 
generally  excluded. 

However,  the  analogue  system,  and  in  particular  the  plate,  allows 
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two  different  characteristic  inodes  of  wave  propagation.  That  is,  for  a 
given  frequency,  the  characteristic  equation,  being  fourth  order  as  is 
the  Orr-Soiranerf eld  equation,  yields  two  branches  corresponding  to  two 
separate  wavelengths  and  phase  velocities.  (The  effect  of  the  foundation 
will  be  to  modify  these  modes  and  most  likely  add  an  additional  charac- 
teristic mode,  or  modes.)  Nevertheless,  in  general  for  a given  excitation, 
or  disturbance,  these  two  modes  are  coupled.  In  other  words,  within  the 
analogue  system  itself  is  an  inherent  mechanism  for  coupling  between 
relatively  long  and  short  wavelength  motions.  For  example,  in  some  cir- 
cumstances, the  mechanism  for  inducing  relatively  short  wavelength 
Tollmien-Schlichting  instabilities  into  a boundary  layer  by  means  of 
relatively  long  wavelength  acoustic  or  structure-borne  disturbances  may 
be  a manifestation  of  this  (linear)  phenomenon.  To  be  precise  the  coupling 
would  require  coincidence  between  the  disturbance  and  one  of  the  stable 
eigenvalues  of  the  Orr-Sommerfeld  equation,  i.e.,  one  that  implies  wave 
attenuation,  which  would  then  be  inherently  coupled  within  the  boundary 
layer  itself  to  an  unstable  mode.  The  extent  to  which  this  represents 
a potential  destabilizing  mechanism  in  connection  with  structure-borne 
noise  in  submarine  hulls  is  beyond  the  scope  of  this  study. 


APPENDIX  A:  THE  RESULTS  OF  A LINEAR  STABILITY  ANALYSIS  INTO  THE  EFFECTS 


OF  A BOUNDARY  WHICH  DISPLAYS  FLEXURAL-LIKE  VIBRATIONS 

1.  Introduction 

14 

Presented  in  this  appendix  are  the  results  of  a study  in  which  the 
inhomogeneous  Orr-Sommerf eld  equation  is  solved  for  the  case  of  a travel- 
ing wave  propagating  along  the  boundary  with  constant  amplitude,  and  real 
wave  number  and  frequency.  In  other  words 

v(x,o)  = -iwC  exp[ik(x  - ct)  1 = — iio^  exp[i(kx  - wt)  1 
o o 

where  v(x,o)  is  the  velocity  imposed,  and  normal  to,  the  boundary  y=0 
and  £ is  displacement  amplitude.  The  results  are  interpreted  in  terms 
of  parameters  relevant  to  the  problem  of  structure-borne  noise  in  sub- 
marine hull  plating. 

2 . Results  of  Analysis 

Using  an  approach  analogous  to  that  used  by  Beniamin* * for  the  homo- 
geneous case  Becker  is  able  to  obtain  explicit  expressions  for  the  wall 
pressure  and  for  the  exchange  of  energy  between  the  boundary  and  the 
adjoining  fluid  for  certain  asymptotic  situations  characterized  by  larqe 
or  small  values  of  the  parameter 

2 

b = ^-  = — 
vk  vw 

For  b small,  these  results  resemble  those  of  the  planar  Couette  flow, 
which  automatically  excludes  Tollmien-Schlichting  instabilities,  and  is 
therefore  uninteresting  as  well  as  unrealistic.  In  contrast,  for  largo 
b,  Tollmien-Schlichting  waves  can  arise. 


(Al) 


(A2) 
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The  inviscid  solution  4^  of  the  Orr— Sommerf eld  equation  yields  the 


time-averaged  power  flow  per  unit  boundary  area,  from  the  wall  into  the 


fluid: 


E-  -Ik  6oI"1+i(o)1 


where  p refers  to  the  density  of  water.  In  Eq.  A3,  f,  is  the  velocity  anpli- 

I I 

tude  of  the  boundary,  and  Im[4^(o))  the  imaginary  component  of  evaluated 
on  the  boundary.  The  energy  flow  associated  with  the  viscid  solution,  being 
small  for  large  b,  need  not  be  discussed.  In  view  of  Eq.  A3,  only  complex 

I 

or  imaginary  values  of  <t>^(o)  give  rise  to  energy  exchange.  This  restriction 
is  satisfied  when  > c,  i.e.  when  a critical  layer  (U(y)  = c)  exists. 


This  requirement  is  also  obtained  from  Benjamin.^  For  kiS  <<  1,  i.e.. 


for  all  situations  of  practical  interest,  it  is  found  that 


4^(0)  = -1Ul) 


4 " 

(l*r  ^ 


where  U and  U are  derivatives  of  U(y)  evaluated  at  the  critical  layer, 
c c 1 


Combining  Eqs.  A3  and  A4,  one  obtains 


4 " 

- - ( U \ U 

TT  2-2  I , “I  c 

= X pc  E u 11-—)  

2 o \ c / . t, 


k6«l 


For  profiles  characterized  by  Uc  < 0,  4^(0)  > 0,  a negative  value  of  E 


results.  This  implies  energy  flow  from  the  fluid  into  the  boundary  which 


is  characteristic  of  all  inviscid  stable  solutions.  For  c > U , 4, (o) 

00  1 


is  real,  and  E = 0. 


Now,  however,  consider  the  case  where  the  parameters  w and  c in 
Eq.  A1  match  the  eigenvalues  of  the  homogeneous  equation  which  implies 


c < U . Under  these  resonant  or  coincident  conditions,  for  the  Blasius 


profile 
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where  is  a Reynolds  number,  Eq.  A2,  corresponding  to  neutral  stability 

of  Tollmien-Schlichting  waves,  say,  w and  k . Since  Eq.  A6  is  positive, 

c c 

this  yields  a negative  value  of  E from  Eq.  A3,  i.e.  energy  flow  from  the 
liquid  to  the  wall.  This  appears  to  be  consistent  with  Benjamin's  solu- 
tion, whereby  Class  A instabilities,  i.e.  Tollmien-Schlichting  waves,  are 
enhanced  if  the  boundary  is  damped. 

Asymptotic  solutions  for  both  kS  <<  1 and  k6  >>  1 indicate  that  energy 
flows  from  the  fluid  into  the  boundary  for  large  Reynolds  numbers.  Becker 
speculates  that  this  is  a property  of  all  velocity  profiles  which  are  stable 
in  the  absence  of  viscosity,  e.g.  the  Blasius  profile.  For  unstable  pro- 
files, i.e.  profiles  displaying  an  "inflexion"  point  (or  a point  where 
dU/dy  becomes  negative)  a different  situation  arises.  Here  energy  can 

I II 

flow  in  either  direction.  For  k6  <<  1 and  kiS  >>  1,  <j>  (o)  < 0 for  U >0, 

1 c 

which  makes  Eq.  A3  positive.  (The  significance  of  this  situation  to  hull 
vibrations  is  discussed  in  Section  IV.  B.) 


APPENDIX  B:  AN  ASYMPTOTIC  LINEAR  STABILITY  ANALYSIS  APPLICABLE  TO  DISTURBANCES 


AND  BOUNDARY  CONDITIONS  ASSOCIATED  WITH  SUBMARINE  HULL  PLATING 

1 . Introduction 

This  appendix  details  the  asymptotic  solution  to  the  Orr-Sommerf eld 
equation  subject  to  the  assumptions  listed  below, 

a.  R >>  1,  i.e.,  Reynolds  numbers  are  large. 

b.  c/Uoo  >:>  i.e.,  the  phase  velocity  of  the  potential  disturbance 

to  the  laminar  boundary  layer  is  much  greater  than  the  free  stream  flow 
velocity. 

Both  assumptions  are  uniformly  valid  for  practical  problems  with  sub- 
marine speeds  and  vibration  spectra  associated  with  propagating  flexural 
wanes.  The  results  of  this  analysis  are  applicable  to  the  physical  situations 
depicted  in  Fig.  5 as  Items  3 and  4,  viz.,  situations  when  the  disturbance  is 
introduced  into  the  boundary  layer  other  than  at  the  boundary. 

2 . Analysis 

For  completeness,  the  Orr-Bommerf eld  equation  is  redefined 

below 


IV  2 " 4 " ? " 

<P  (y)  -2a  <t>  (y)+a  <t>(y)  = iaR{  (U-c)  [<|>  (y)  -a'  <f>  (y)  ] -U  4>(y)}  (Bl) 

<p(y)  is  defined  by  the  perturbation  stream  function  <I>(x,y)  which  is 
assumed  to  be  of  the  form 

4>(x,y)  = 4>  (y)  exp[ia(x  - ct)  1 (B2) 

x and  y are  cartesian  coordinates  parallel  and  normal  to  the  direction 
of  flow.  R is  the  Reynolds  number.  All  equations  are  in  non-dimensional 
form  whereby  lengths  are  normalized  to  the  boundary  layer  thickness  6 ; 

velocities  to  the  free  stream  velocity  U , time  by  fi/U  and  pressure  by 

2 un  °° 

pU  where  p is  the  fluid  density. 

00  • 

The  stability  problem  is  thus  represented  by  the  solution  to  Eq.  Bl 
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subject  to  appropriate  boundary  conditions. 

Consider  the  case  of  parallel  flow  over  a flat  boundary  of  infinite 
extent.  The  appropriate  conditions  at  the  outer  edge  of  the  boundary  layer 
can  be  expressed  by 


>(“>)  = <j)  (<x>)  = o 


(B3) 


At  the  boundary,  continuity  of  velocity  tangential  to  the  boundary  dictates 
that9 


(1  - c/Z  ) [<f>  (0)  + (U  (0)/c)<(>(0)] 


+ i(ZtaR)_1  [<*>'"  (0)  - aV(0)]  = 0 


(B4 ) 


where  Z^  is  the  tangential  impedance  of  the  boundary. 

for  the  fourth  boundary  condition  continuity  of  normal  velocity  is 
prescribed  in  terms  of  the  normal  impedance  of  the  boundary 


p(0)/  [-ia<J>  (0)  ] = -Z 


(B5) 


where  p(0)  is  the  perturbation  pressure  at  the  boundary. 

Boundary  conditions  (B4)  and  (B5)  pertain  to  the  modifying  effect  of 
a finite  impedance  hull  on  disturbances  in  the  laminar  boundary  layer.  The 

x-momentum  equation  can  be  conveniently  used  to  express  p(0)  in  terms  of 

★ 

the  stream  function  . 


p (0)  = — i [ 4>  (0)  - a <f>  (0)  ] / (aR)  + c<t>  (0)  + U (0)<j>(0) 


(B6) 


Consider  the  behavior  of  Eq.  B1  for  large  values  of  R,  i.e.  R >>  1. 
In  this  range,  Eq.  B1  can  be  satisfied  via  two  mechanisms: 


* 13 

This  expression  differs  from  that  used  by  Benjamin  which  was 

derived  using  the  y-momentum  equation  and  is  of  integral  form. 
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a.  <p  exhibits  large  gradients  thus  making  the  left  hand  side  of 
Eq.  B1  comparable  to  the  right  hand  side  or 

b.  | (U  - c)  (<p"  - a2<j>)  ~ <<  1. 

Mechanism  (a)  is  a cause  of  considerable  difficulty  in  developing  numerical 
integration  solutions  to  Eq.  Bl. 

Eq.  Bl  with  R » 1 is  of  the  form  of  a singular  perturbation  problem 
in  that  the  highest  ordered  derivative  term  is  multiplied  by  a small  number, 
namely  eSR  ^ . This  can  result  in  large  gradients  being  exhibited  by  the 
solution  in  order  to  satisfy  imposed  boundary  constraints.  The  classical 
approach  to  such  problems  is  to  define  an  "innei"  solution  that  is  valid 

within  the  vicinity  of  those  large  gradients,  and  to  "match"  this  solution 

. 19 

with  an  "outer"  solution  that  is  valid  beyond  this  region.  The  particular 

problem  under  consideration  is  somewhat  degenerate  in  that  the  first  order 

outer  solution  which  is  obtained  by  setting  R = ® (or  f = 0)  in  Eq.  Bl, 
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is  identically  zero  assuming  the  quantity  ca  <<  1.  This  means  that  the 
inner  solution  must  be  uniformly  valid  across  the  boundary  layer. 

In  order  to  obtain  the  appropriate  inner  solution,  Eq.  Bl  is  subjected 
to  a change  in  the  independent  variable  in  the  form  n = y/en.  The  resulting 
equation  is  then  solved  via  a power  series  in  e.  We  will  only  deal  here 
with  the  first  order  approximation. 

If  the  term  (U  - c)  changes  sign  within  the  domain  of  the  boundary 
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layer  the  appropriate  transformation  is  given  by  n = 1/3.  This  yields 
a solution  in  terms  of  Bessel  functions  of  order  1/3,  and  in  turn  the 
familiar  Tietjens  function.  However,  for  the  case  at  hand,  where 
| U/c | <<  1 throughout  the  boundary  layer,  a more  convenient  transforma- 
tion results  when  n = 1/2.  Thus  the  "stretch"  variable  n is  defined  by 

g = y/e ^ (B7) 

Substituting  Eq.  B7  into  Eq.  Bl  and  allowing  f=>0  yields  the  equation 

*nnnn(n)  + ic%n(n)  = 0 (c2)  5 0 m) 

The  general  solution  to  Eq.  B8  is  given  by 


♦ = Asin [ (1  + i) an]  + Bcos[(l  + i)an]  + C + Dn  (B9) 

with  a = (c/2)  . 

However,  Eq.  B3  represents  a severe  constraint  on  Eq.  B9  and  yields 
a non-trivial  solution  only  if 

C = D = 0 (B10) 

and  either 

B = "iA  (Blla) 

or 

Re  [a]  = -Im [a]  (Bllb) 

Condition  Bllb  implies  that  c is  imaginary  and  negative  and  this,  in  turn, 
implies  a stable  solution. 

Thus,  the  problem  is  reduced  to  the  task  of  determining  if  a solution 
in  the  form 

<t>  = A[sin(l  + i)  an  - icos(l  + i)ag]  (B12) 

is  compatible  with  the  boundary  conditions  represented  by  Eqs.  B4  and  B5  and, 
if  compatible,  whether  an  unstable  mode  is  possible.  The  fact  that  Eq.  B12 
exhibits  only  one  free  constant.  A,  requires  a compatibility  condition  on 
the  boundary  conditions  themselves. 

A solution  exists  only  if 

zt  = ■*„♦«»/(  (o' <0)/c>*«»  + <f>n(0)/e]  (B13) 

'0  c -+  0 

Now,  subjecting  Eq.  B6  to  the  coordinate  transformation  given  by  Eq.  B7 

33 


yields 


ep(0)  = -i  [<t>  (0)  - a2e2<(>  (0))  + c<(>  (0)  + Fu ’ (0)  4>  CO)  (B14) 

nn  n n o 

Using  Eq.  B8  and  dropping  c terms  yields  the  result 

p(0)  = u'  (0)  <)>  (0)  (B15) 

which  when  substituted  into  Eq.  B5  gives 

Z = -iU  (0) /a  (B16) 

n 

If  Eq.  B16  is  to  be  interpreted  in  dimensional  form,  then  the  right  hand 
side  of  the  equation  must  be  multiplied  by  p. 

3.  Discussion  of  Results 

Eqs.  B13  and  B16  give  the  values  of  the  boundary  impedances  that  are 
compatible  with  a neutrally  stable  condition.  It  should  be  noted  that  for  a 
disturbance  traveling  in  the  same  direction  as  the  flow,  i.e.  (Re(a)  > 0), 
the  required  normal  impedance  is  inertial.  However,  for  a disturbance 
traveling  opposite  to  the  flow,  the  impedance  is  stiffness-like.  The  effect 
of  energy  dissipation  within  the  boundary  can  also  be  analyzed  from  Eq.  B16. 
Boundary  dissipation  would  correspond  to  a positive  real  component  of  Z^. 
However,  from  Eq.  B16  this  would  require  a negative  imaginary  component  of 
a which  implies  an  unstable  condition  when  Re (a)  > 0 and  a stable  condition 
when  Re(oi)  < 0 (Eq.  B2) . Therefore  the  effect  of  boundary  dissipation  is 
destabilizing  for  potential  disturbances  traveling  in  the  direction  of 
flow  and  stabilizing  when  the  disturbance  travels  against  the  flow.  More 
importantly  condition  B13  which  is  required  for  an  instability  is  not 
satisfied  by  present-day  submarine  hull  construction. 
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Enveloped  phase  velocities  vs.  frequency  plot  for  freely 


Figure  4 Spatial  Amplification  Rates  for  the  Blasius  Boundary  Layer 

over  a Rigid  Surface  k x R = (cr/0j(a  y/Uj],  Reproduced  from  Reference 


Potential  Destabilizing  Effects  of  Structure-Borne  Noise  on  Hypothetical  Laminar  Flow  Submersible 


FIG.  6 TIMOSHENKO -MINDLIN  PLATE  RESTING  ON  A LOCALLY 
REACTING  FOUNDATION 
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